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ABSTRACT 

Elastic scattering on nucleons, vN — > Nu, is the dominant supernova 
(SN) opacity source for \i and r neutrinos. The dominant energy- and 
number-changing processes were thought to be ve~ — > e~v and vv <-> e + e~ 
until Suzuki (1993) showed that the bremsstrahlung process vvNN <->• NN 
was actually more important. We find that for energy exchange, the related 
"inelastic scattering process" vNN NNv is even more effective by about 
a factor of 10. A simple estimate implies that the and v T spectra emitted 
during the Kelvin-Helmholtz cooling phase are much closer to that of v e than 
had been thought previously. To facilitate a numerical study of the spectra 
formation we derive a scattering kernel which governs both bremsstrahlung and 
inelastic scattering and give an analytic approximation formula. We consider 
only neutron-neutron interactions, we use a one-pion exchange potential in 
Born approximation, nonrelativistic neutrons, and the long- wavelength limit, 
simplifications which appear justified for the surface layers of a SN core. We 
include the pion mass in the potential and we allow for an arbitrary degree of 
neutron degeneracy. Our treatment does not include the neutron-proton process 
and does not include nucleon-nucleon correlations. Our perturbative approach 
applies only to the SN surface layers, i.e. to densities below about 10 14 gcm~ 3 . 
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1. Introduction 

A quantitatively accurate prediction of the fluxes and spectra of supernova (SN) 
neutrinos is required for a theoretical understanding of an array of fascinating phenomena, 
including the explosion mechanism itself, the cause of neutron star natal kicks, various 
aspects of SN nucleosynthesis, and the interpretation of the neutrino signal from SN 1987A 
and future galalactic SNe, with or without the assumption of neutrino masses and mixings. 
However, at the present time one is far away from this goal because of numerical limitations 
and because of significant shortcomings in the calculation of the relevant microphysics, 
notably the equation of state and the neutrino opacities. 

We presently study the opacity contribution of nucleon bremsstrahlung AW — > NNvv 
and its inverse uuNN — > AW, processes which have been ignored in all numerical SN 
studies except for the proto-neutron star cooling calculations of Suzuki (1991, 1993) who 
found a significant modification of the neutrino fluxes and spectra. We also study the 
inelastic scattering process uNN — > NNv which is obtained when crossing a final-state 
bremsstrahlung neutrino into the initial state. Janka et al. (1996) had stressed its apparent 
importance for the neutrino spectra formation, but without making a connection with 
Suzuki's work. We will again motivate the importance of these processes and provide a 
scattering kernel which allows for a practical implementation in a numerical SN code. 
While this is a modest aspiration relative to the large number of open questions regarding 
SN neutrino opacities, we hope that our work may nevertheless prove useful for future 
numerical studies. 

In order to appreciate the importance of nucleon bremsstrahlung and related reactions 
recall that for u e and u e the dominant opacity contribution derives from charged-current 
processes {j3 processes) which involve electrons and positrons. In the diffusion regime they 
keep v e and v e essentially in local thermal equilibrium. We are mostly concerned, however, 
with the transport of u^, u T , and their antiparticles, to which we will collectively refer 
as v^. The dominant opacity contribution is the neutral-current scattering on nucleons 
v^N — > Nu^, a process which does not equilibrate the neutrino number density and which 
is ineffective at modifying the spectrum because the neutrino energies are low relative 
to the mass of the nonrelativistic nucleons, especially in the outer layers of a SN core. 
Therefore, processes such as v^e~ — > e _ z/ M and v^v^ <-> e + e _ , which are subdominant 
with regard to the total opacity, are nevertheless important for the equilibration of the 
neutrino number density and spectra. It turns out that nucleon bremsstrahlung is far 
more effective than pair annihilation at equilibrating the neutrino number density, and it 
is of comparable importance to f M e~ — > e~z/ M at equilibrating the spectra. Moreover, the 
inelastic scattering process vNN NNv is in turn far more effective at exchanging energy 
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than bremsstrahlung and is thus the dominant process for the neutrino spectra formation. 
Including these processes in a proto-neutron star cooling calculation would make the 
spectrum far more similar to that of u e than had been thought. 

While nucleon bremsstrahlung and related processes are conceptually simple, a reliable 
calculation is nevertheless nontrivial. Difficulties include the nuclear matrix element, i.e. 
the appropriate nucleon-nucleon interaction potential, the intermediate degree of nucleon 
degeneracy, NN correlations, and the role of multiple-scattering effects. We will not be 
able to resolve all of these issues. For example, we will completely ignore NN correlations 
which may be quite important in some regions of the SN core where nuclei may not even be 
completely dissociated. In regions where nuclei exist, nuclear bound-bound or bound- free 
transitions involving neutral-current neutrino reactions may be important. The main 
advance of our calculation is the treatment of intermediate degrees of nucleon degeneracy 
in free-free transitions, the inclusion of the pion mass in the nucleon interaction potential, 
and the inclusion of multiple-scattering effects. Our scattering kernel is then self-consistent 
in the sense that it allows for the simultaneous treatment of bremsstrahlung and inelastic 
scattering while producing the correct total uN scattering cross section. 

We focus on the neutrino spectra formation which takes place in the surface layers of a 
proto neutron star where number-equilibrating, energy-equilibrating, and finally scattering 
processes freeze out at different radii so that the different neutrino species are emitted with 
different fluxes and different spectral temperatures, but approximately equal luminosities. 
We expect that bremsstrahlung and inelastic scattering make the spectra and fluxes more 
similar than had been thought previously. 

Nucleon-nucleon interactions not only lead to bremsstrahlung and the inelasticity of 
neutrino-nucleon scattering, but also to a reduction of the overall uN scattering cross 
section (Raffelt & Seckel 1995; Sawyer 1995; Raffelt, Seckel & Sigl 1996). This cross-section 
reduction is an inevitable consequence of the fact that nucleon spins interact so that one 
cannot consistently include bremsstrahlung and the inelasticity of uN scattering and yet 
ignore this cross-section reduction. Unfortunately, for densities above, say, 10 14 gcm~ 3 a 
controlled calculation of the scattering kernel and thus of the neutrino transport coefficients 
is currently out of reach. The only hint as to the possible magnitude of the cross-section 
reduction in the inner parts of a SN core derives from the signal duration of SN 1987A 
(Keil, Janka & Raffelt 1995) and the /-sum rule of the spin-density structure function (Sigl 
1996) which imply that a naive extrapolation of the perturbative spin-density structure 
function into the high-density regime would significantly overestimate the cross-section 
suppression. Recent 2-dimensional hydrodynamic SN cooling calculations indicate that 
actually convection, not neutrino radiative transport, may be the dominant mode of energy 
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transfer from the inner proto-neutron star to the neutrino sphere (Keil, Janka & Miiller 
1996), rendering the high-density opacities less critical for the overall neutrino luminosity. 
In any case, the problem of the neutrino spectra formation in the surface layers is rather 
disjoint from the problem of calculating the overall neutrino luminosity. 

We begin our discussion in Sec. 2 with a comparison of the different neutrino opacity 
contributions and thus motivate the importance of nucleon bremsstrahlung and related 
processes. In Sec. 3 we formulate the bremsstrahlung rate in terms of the dynamical 
spin-density structure function (the scattering kernel) and thus reduce the problem to the 
calculation of a dimensionless function which includes all of the nuclear and many-body 
complications. We also discuss the underlying nuclear matrix element of the bremsstrahlung 
process. In Sec. 4 we study the nucleon phase space and provide the appropriate analytic 
approximation formulae. In Sec. 5 we summarize and discuss our results. 



2. Comparison of Opacity Sources 

2.1. Mean Free Path 

In a SN core the main opacity source for \i and r neutrinos is neutral-current scattering 
on nucleons. However, the large nucleon mass relative to typical neutrino energies renders 
this process ineffective at equilibrating the neutrino spectra and it certainly cannot modify 
the neutrino number density. We begin with a comparison of the relative importance of 
those processes which are subdominant with regard to the total opacity, yet allow for an 
effective modification of the spectrum and density. The simplest case is v^v^ — > v e v e . 
Using Maxwell-Boltzmann distributions at temperature T for all neutrino species, we find 
the thermally averaged absorption rate 

r = 4 G ^ 5 , (i) 

7T° 

which sets a natural scale for all other processes. We use natural units with h — c — k& — 1. 
The rate for z/^z/ e — > f M z/ e together with that for z/^^ — > z/^z/ e is 4r . For u^i?^ —>■ e + e~ it 
is T (C£ >e + C\ e )if e e-^ jYl while it is T (C* e + C% e )^ 2 e for -> z/ At e ± , where rj e > 1 
is the electron degeneracy parameter. The weak coupling constants are Ca,s = — \ and 
Cv, e = — \ + 2sin 2 w for and v T . Neutrino coalescence vv — > plasmon is negligible. 
Finally, there is inverse bremsstrahlung on nucleons v^v^NN — * NN\ the spectrally 
averaged absorption rate is given by Eq. ( |25|) in Sec. 3.4. 

In order to compare these different processes for realistic conditions we use the 
numerical SN model S2BH_0 of Keil, Janka & Raffelt (1995) which represents a SN core 
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1 s after collapse; the radial profiles of various physical parameters are shown in Fig. 1. In 
the outer parts of the star, which is where the neutrino spectra are formed, the effective 
nucleon degeneracy parameter 77* = p 2 F /2m^T varies between 2 and 4, n e between 1 and 8, 
and the range of relevant temperatures is taken to be T = 5-10 MeV, leading to 



r = r n x <^ 
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(2) 



z/ M e + plus v^- 



We conclude that for typical SN conditions bremsstrahlung is by far the most important 
number-changing reaction. Even for energy exchange it is more important than elastic 
scattering on electrons and positrons. 



2.2. Energy Transfer 

The processes considered in the previous section are very effective at transferring 
energy in the sense that in a given interaction the neutrino essentially loses all memory of 
its previous energy. This is not the case for the "inelastic scattering process" vnn — > nnv 
which is the bremsstrahlung process with a final-state neutrino crossed into the initial 
state. While we use the term "inelastic scattering process" we stress that it is not logically 
distinct from the "elastic channel" vn — > nv. Both are described by one and the same 
scattering kernel which has a finite width as a function of the energy transfer uj due to the 
nucleon-nucleon interaction. Therefore, in each collision the final-state neutrino energy is 
smeared out by a small amount which is given by the width of the scattering kernel. 

To quantify the efficiency of energy transfer of this process relative to those of the 
previous section it is useful to consider a fluid of neutrinos at a temperature T v which is 
slightly different than the temperature T of the neutron bath. According to Eqs. fl3"T| ) and 
([34]) we find for the rate of energy transfer between these fluids 

AQ _ 3C 2 A G 2 F n B T 5 T a f 1 for nn <-> nnvv, 
AT ~~ tT 3 2ttT X I 20 for ran «-> nnv, 

where the spin-fluctuation rate r CT is given in Eq. (]20|) . For the bremsstrahlung case we 
have included a factor 1/2 relative to Eq. (|31|) ; we only count the energy transferred to the 
neutrinos. These results receive numerical corrections of order unity for realistic conditions, 
reducing the relative importance of scattering to perhaps a factor of 10. 
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For completeness we have calculated the same quantity for neutrino scattering on 
degenerate electrons, 

AQ 2^ = ^ Gl{Cl e + Cfarfc*, (4) 
where again Cy, e ~ and CA, e = — \ for v^. Relative to inelastic nucleon scattering this is 

AQ^ e -^ e -^ = 1 ^ x io _ 3 r£ T y 2 ^ (5) 

Vf 1 nn*-*nnv fl V* 

where T 10 = T/10 MeV. In the outer part of our SN model we have r\ e < 8, r)* = 2-4, and 
Tio = 0.5-1 so that the ratio is always smaller than 0.01. Electrons are negligible for the 
energy transfer, in complete keeping with the conclusions of Janka et al. (1996). 

In the calculation of bremsstrahlung and related processes we assume nonrelativistic 
nucleons and use the long-wavelength approximation. As a consequence neutrinos cannot 
transfer energy to the nucleons if the NN interaction is switched off. From Janka (1991) 
and Tubbs (1979) and using a Maxwell-Boltzmann distribution for the neutrinos we find 

AQrecoii = (C v + 5C 2 A )G 2 F n B T 5 240 T 
AT 7T 3 m N ' U 

Naturally, the vector-current interaction appears here in contrast to bremsstrahlung-related 
effects. Ignoring Cy we find 



^Qvnn^nnv _ f^N^a _ 2 2y 27T 3/2 _ 91 3/2 /y\ 



where we have used Eq. (|20j) for the spin-fluctuation rate. The vector current contribution 
to recoils reduces this ratio a bit, and pion-mass effects reduce the inelastic term by another 

3 /2 

small amount so that it may be more realistic to take something like 0.5 77* . In the SN core 
model of Fig. 1 we always have 77* = 2-4 so that recoil effects are always of roughly equal 
importance to inelastic scattering. Therefore, the transfer of energy is even more effective 
than what is accounted for by our scattering kernel which ignores nucleon recoils. 



2.3. Energy Sphere 

We may next attempt to estimate the change in the spectral temperature of the 
emitted neutrino fluxes which is brought about by including bremsstrahlung and related 
processes. Neutrinos stream off freely from the "transport sphere" which is the radius 
where scattering on nucleons is no longer effective at trapping them. Deeper inside is the 
"energy sphere" which is the radius where energy-exchanging processes such as scattering 
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on electrons become ineffective (Burrows & Mazurek 1982, 1983). We presently estimate the 
temperature of the energy sphere for the SN model of Fig. 1 with or without the inclusion 
of bremsstrahlung and inelastic scattering. We take the change in this temperature as 
an indication of the changed spectral temperature of the emitted neutrinos. We stress, 
of course, that the temperature at the energy sphere is not identical with the spectral 
temperature of the emitted neutrino flux which is modified by the fact that lower energy 
neutrinos diffuse more effectively through the layer between the energy and the transport 
spheres. Still, we think that the following procedure gives one a sense of the quantitative 
importance of bremsstrahlung and its sister reactions. 

In order to calculate the location of the energy sphere we introduce the effective mean 
free path for energy exchange (Shapiro & Teukolsky 1983) 



Here, A to t is the mean free path for a neutrino of a given energy while A e is the mean free 
path against energy-changing reactions. If energy-changing reactions dominate the total 
opacity we have A to t — A e and thus A c g- ~ A c . However, in a "scattering atmosphere" as in 
a SN core where A to t — A scat tcr the trapping of neutrinos by collisions gives them a larger 
chance to lose energy by some other process, leading to an increased value of A^ 1 relative 
to A" 1 . We define the radius R e of the energy sphere by 



We have performed this calculation for the SN model of Fig. 1 first for muon 
neutrinos by including the neutral-current scattering on protons and neutrons as the 
dominant scattering process. Nucleon degeneracy effects were taken into account. As 
an energy-exchange process we first use electron scattering, again including degeneracy 
effects. We then find that the temperature of the energy sphere varies with neutrino 
energy as shown by the solid line in Fig. 2. Next we use inverse bremsstrahlung as the 
only energy-changing reaction according to the numerical prescription developed in Sec. 4, 
leading to the short-dashed line. We see that the two processes are of roughly equal 
importance, with electron scattering dominating for relatively large neutrino energies. We 
have repeated the same exercise for z/ e , except that we include the charged-current reaction 
u e p — > ne + as an energy- changing reaction, leading to the long-dashed curve. Evidently this 
is the most important reaction for u e and thus dominates the spectrum formation. 




(8) 




(9) 



i.e. by the radius where the effective optical depth is |. 



To extract one fixed energy-sphere radius we next take neutrinos to be in thermal 
equilibrium up to their energy sphere so that their distribution is characterized by the local 
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medium temperature, and that further out they are characterized by their energy-sphere 
temperature. All reaction rates are averaged with this spectral distribution so that 

(Ki) = ^(KoDiK 1 )- (io) 

The energy sphere is defined as in Eq. @ except using (X~g) instead of A^ 1 . In order to 
locate the neutrino sphere we must now perform an iteration, assuming first some estimate 
for the energy-sphere temperature to determine the location where the effective optical 
depth is |, then use the temperature there as the next approximation, and so forth until 
convergence (Keil & Janka 1995). We stress that our benchmark model was calculated with 
equilibrium neutrino transport everywhere so that it becomes actually more self-consistent, 
not less so, when the neutrinos are kept in equilibrium out to larger radii. 

For the SN core model of Fig. 1 we find the energy-sphere temperatures and densities 
shown in Table 1. We have included the bremsstrahlung process according to the 
prescription developed in Sees. 3 and 4, and then multiplied its rate by a factor /brans- 
For /brems = we find the energy sphere without bremsstrahlung while /brems = 1 includes 
the full effect. We also show the redshifted energy-sphere temperatures for an observer at 
infinity and the total optical depth at the location of the energy sphere, again for a spectrally 
averaged neutrino mean free path. We are thus led to conclude that bremsstrahlung alone 
reduces the spectral temperature by about 1 MeV, while that of v e remains largely 
unchanged. The change of the spectrum agrees with Suzuki's (1991, 1993) results. 

In order to estimate the impact of inelastic scattering we recall that it is by a factor of 
10 more important than bremsstrahlung, give or take a factor of 2. Therefore, we have also 
calculated the energy sphere for /brems = 5, 10, and 20, which probably mimics the effect 
of inelastic scattering and its uncertainty. The v e temperature is still only mildly affected, 
suggesting that (3 processes are still very important for this species. The z/ M temperature is 
significantly lowered; it is very close to that of v e . 

Of course, including NN interactions consistently will probably accelerate the cooling 
process due to the reduced overall neutrino scattering cross section in the deep interior of 
the SN core. This will likely heat the neutrino sphere so that the spectral temperatures may 
actually increase. A similar effect will be caused by convection if it is a generic phenomenon 
for the proto-neutron star evolution. We believe, however, that our schematic treatment 
gives us a reasonable estimate of the differential temperature change, i.e. the and v e 
temperatures will be much closer than is often assumed. 

A quantitative assessment of the spectral modifications requires detailed numerical 
simulations which we are in no position to perform. In the following, however, we provide a 
scattering kernel which is needed for such an investigation. 
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3. Bremsstrahlung and Related Processes 

3.1. Scattering Kernel 

For a numerical study the bremsstrahlung process, its inverse, and inelastic scattering 
are best formulated in terms of a "scattering kernel" which contains all of the properties of 
the nuclear medium, but which excludes the neutrino phase-space integration. We begin 
with the usual weak interaction Hamiltonian density^] 

G F - 

'Hint = ~j= i'Nl^Cy - CaJbJ^N^uJ^ - lb)^u (11) 

where Gf is Fermi's constant, ?/V the nucleon Dirac field for either protons or neutrons, 
and ip u the neutrino field. The neutral-current vector coupling constant is \ — 2 sin 2 0w ~ 
for protons and — | for neutrons. The axial-current coupling is often taken to be ±1.26/2 
for protons and neutrons, respectively, but the strange-quark contribution to the nucleon 
spin causes certain deviations from this simple picture (Raffelt & Seckel 1995). 

We are concerned with conditions where the temperature is below 10 MeV so that a 
typical nucleon velocity is around 0.2 c. Therefore, it is reasonable to use the nonrelativistic 
limit which implies that only the axial vector current contributes to bremsstrahlung (Friman 
& Maxwell 1979; Raffelt & Seckel 1995). The squared matrix element for the emission of a 
neutrino pair, N\ + N2 — > N% + N4 + v + P, may be written as 

£|.M| 2 =(^)ViV- (12) 

spins \ V / 

where M^ v and N^ v stand for the nucleonic and neutrino parts, respectively. From the 
discussion in Raffelt and Seckel (1995) it follows that in an isotropic medium and for 
nonrelativistic nucleons it is enough to consider the spatial trace M = | M\ which defines 
the "reduced squared matrix element." Equation (|TJ) may thus be written as 

]T \M\ 2 -> (2^L\ M8aW3-cos0), (13) 

spins \ V & / 

where u\ t 2 are the neutrino energies and 9 is the angle between their momenta. With these 
simplifications and in the long-wavelength limit the scattering kernel is 

n B J H[ 2mjv(27r) 

x (2n) 4 5(E 1 + E 2 -E 3 -E 4 + u) <5 3 ( Pl + p 2 - p 3 - p 4 ), (14) 



1 In previous papers (e.g. Janka et al. 1996) it had been written (Gf/2v2) ... so that, for example, Ca 
was ±1.26 rather than ±1.26/2 which we use here. 
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where n# is the baryon density. Further, /j is the occupation number of nucleon i with 
energy Ei and momentum pj while u = —(uji + uj 2 ) for pair emission and oo = iO\ + cj 2 for 
pair absorption, i.e. u; is the transfer of energy to the nuclear medium. For neutron-neutron 
or proton-proton bremsstrahlung we need to include a statistics factor 1/4 to compensate 
for initial- and final-state double counting. The superscript (1) signifies that in this form 
the scattering kernel is a lowest-order perturbative result, derived from a bremsstrahlung 
calculation. The subscript a indicates that in the framework of linear-response theory it is 
the nucleon spin-density autocorrelation function (e.g. Janka, Keil, Raffelt & Seckel 1996). 
Both the lowest order S^(u) and the full spin-density autocorrelation function S (7 (lu) obey 
detailed balance S a (— uo) = S a (uo) e^l T . 

In our treatment we always assume that the neutron spins evolve independently of each 
other, i.e. we ignore possible spin-spin correlations. In this case the full scattering kernel 
must obey the normalization requirement (Raffelt & Strobel 1997) 



where f p is the occupation number of a neutron with momentum p and the factor 2 
represents the two spin orientations. In the nondegenerate case where we may neglect the 
Pauli blocking factor (1 — f p ) the normalization is unity. 



In order to calculate the scattering kernel we need to know the nuclear matrix 
element. The bremsstrahlung emission of neutrino pairs or axions arises from nucleon spin 
fluctuations in collisions so that one needs a spin-dependent nucleon-nucleon potential. 
The most general velocity-independent interaction has a scalar (spin-independent) part, a 
central part proportional to <T\ • cr 2 , and a tensor part proportional to 3f • <j\ r ■ <x 2 — <J\ • <r 2 
where <x 12 are the spin operators of the two nucleons (Blatt & Weisskopf 1979, Ericson & 
Weise 1988). These interactions have the important property that they conserve the square 
of the total nucleon spin, i.e. (cti + cr 2 ) 2 is a constant of the motion. In addition, the scalar 
and central parts conserve the total spin (er 1 + cr 2 ). 

We first consider bremsstrahlung from the collision of identical nucleons (proton-proton 
or neutron- neutron collisions). The emission of neutrino pairs or axions is induced by the 
total spin operator (<7i + cr 2 ) which evolves nontrivially only due to the tensor interaction. 
Because the physical cause of the tensor interaction is primarily one-pion exchange (OPE) 




(15) 



3.2. 



Nuclear Matrix Element 



3.2.1. Identical Nucleons 
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one expects that an OPE ansatz captures the dominant aspect of the bremsstrahlung 
process (Friman & Maxwell 1979). In addition, the spin operator has nonvanishing matrix 
elements only between triplet states. Because the total wavefunction must be antisymmetric 
the orbital angular momentum of the relevant states must be odd. Therefore, s-waves do 
not contribute so that the least well-known short-distance part of the nucleon interaction 
potential does not affect the nucleon motion, further justifying the OPE ansatz. 

For identical nucleons the OPE squared matrix element in its "reduced form" is found 
to be (Raffelt & Seckel 1995) 



64(4vr) 2 al 
3 uj 2 



q 2 \ , ( q* \ , q 2 q*-3(q-q* 



q 2 + m 2 J yq 2 + m 2 / (q 2 + m 2 )(q 2 + m 2 ) 



(16) 



where u = E± + E% — E% — E\ is the energy transfer to the nucleons, q = P2 — P4 is the 
momentum transfer between the nucleons, and q* = P2 — P3 the momentum transfer for the 
exchange amplitude. Further, a n = (f2mi^/m n ) 2 /4:'ir pa 15 with / « 1 is the pion-nucleon 
"fine-structure constant." 

At a temperature T a typical nucleon momentum is (3mjyT) 1//2 so that a typical 
momentum exchange in a collision is of a similar magnitude. At T = 10 MeV this is about 
170 MeV, only slightly larger than the pion mass of 135 MeV. First, this implies that the 
pion mass cannot be ignored in the denominators in Eq. ([T6|) . Secondly, it means that the 
typical potential region probed is not much smaller than m~ l so that, again, the OPE 
potential should be a reasonable approximation. 

Two-pion exchange effects become important at distances below 2fm ~ 1.5 m" 1 
(Ericson & Weise 1988). Their impact on the bremsstrahlung process can be estimated by 
mimicking the two-pion exchange contribution by one- p- meson exchange where the mass 
of this effective particle is taken to be m p pa 600 MeV. A typical term in Eq. fllq) is then 
modified to (Ericson & Mathiot 1989) 

' q w q ,! (it) 



y q z + m%. ) \ q + q 1 + m 

with C p = 1.67. Taking q 2 to be 3Tm^ with T = 10 MeV, i.e. taking |q| to be around 
170 MeV yields a 35% reduction of the squared matrix element. This estimate quantifies 
the error one is likely to make by using a simple OPE potential. 

Even if the OPE potential is taken to be appropriate, this alone does not guarantee that 
the Born approximation is justified which was used to calculate Eq. fll6|). Ericson & Weise 
(1988) compare the p-wave scattering volumes of an OPE Born calculation with an iterated 
OPE calculation and with a full calculation using the Paris nucleon-nucleon potential. The 
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OPE Born approximation, again, yields a certain overestimate so that for our conditions of 
interest the correct result may be smaller by as much as 50%. These uncertainties set the 
scale of precision that one may hope to achieve with our simple approach. 

The OPE approximation may be tested by comparing the calculated rate for the 
pionic bremsstrahlung process pp — > ppn° with experimental data (Choi, Kang & Kim 
1989; Turner, Kang & Steigman 1989). The pion coupling is of derivative nature so that 
only the axial-vector current contributes, rendering this process a good proxy for the 
bremsstrahlung emission of neutrino pairs. However, because of the pion mass threshold 
the kinematical regime probed is always at higher (but not very much higher) energies than 
what is appropriate for a SN core, and one needs to include relativistic corrections which 
are no longer negligible. The agreement between an OPE calculation and the data found 
by Turner, Kang & Steigman (1989) is far better than one would have expected according 
to our above caveats. After all, one now probes the NN interaction potential at distances 
where two-pion exchange and other corrections surely must be important. This example 
illustrates the well-known fact that the OPE potential in Born approximation often yields 
far better results for spin-dependent processes than one is entitled to expect. 



3.2.2. Proton-Neutron Scattering 

The situation is far more complicated for proton-neutron scattering. Because the 
neutrino neutral-current coupling to protons and neutrons is nearly equal but of opposite 
sign, the pair emission is now induced essentially by the operator (j\ — er 2 which is not 
conserved by the central part of the interaction potential. This operator connects triplet 
and singlet states, i.e. the selection rules of the neutron-neutron system do not apply. 
Further, the two-nucleon wave function can be antisymmetric in the isospin variables so that 
orbital s-waves allow for both singlet and triplet states. Put another way, the short-distance 
behavior of the potential does matter in this context. By the same token, the comparison 
of Turner, Kang & Steigman (1989) with experimental data is not directly relevant. 

Recently Sigl (1997) studied numerically the behavior of S a (uj) for this case, using 
a phenomenological potential adapted to fit low-energy scattering data and deuteron 
properties. He found huge differences between the exact S a (u) and the one derived in Born 
approximation. Put another way, because s-wave scattering dominates in this problem, 
the phase shifts are large and the Born approximation is not justified. For the chosen 
temperature and density, processes of the form p + n — > d + u + u play a dominant role while 
below the deuteron binding energy (\uj\ < 2.2 MeV) only free-free transitions contribute 
to S a (uj). In this regime Sigl finds that his numerical S a (cu) and the one derived from the 
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OPE potential in Born approximation agree surprisingly well — there must be compensating 
effects between a potential and an approximation which are separately unjustified. 

Thus one may probably estimate the relative significance of the pn process by using 
the OPE Born expression. One finds for the squared matrix element if one takes \Ca\ to be 
the same for protons and neutrons (Raffelt & Seckel 1995) 

"(-ifr)' + 2 ( -AT + ^ . ff ~^ 2 fJ (is) 

with all coefficients the same as in Eq. ([16]). In the limit of a vanishing pion mass the 
expression in square brackets of Eq. ([16]) reduces to 3 — (3 where (3 is a phase-space average 
of 3(q ■ q*) 2 /q 2 q^ which for nondegenerate neutrons is found to be (3 ~ 1.31. In the present 
case the expression in square brackets reduces to 5 — 2/9/3. If we take for (3 the same value 
as before we find that the pn squared matrix element is about 2.5 times as large as the 
nn one. For the np process the bremsstrahlung rate is proportional to the product of the 
densities n p n n , i.e. to n 2 B /4 for an equal mix of protons and neutrons. For pure neutrons 
it is also proportional to n 2 B /4 where the factor 1/4 now derives from the phase-space 
reduction for identical particles. Put another way, the relative importance of the two cases 
is indeed well estimated by a comparison of the squared matrix elements. 

For our conditions of interest protons are relatively rare so that ignoring the pn process 
leads to an underestimate of the bremsstrahlung rate which is probably not larger than a 
few tens of percent. Therefore, the error made by considering a medium of neutrons alone 
partially compensates the error made by using the OPE potential in Born approximation. 
We are thus led to believe that the compound error of our approximations does not exceed 
a few tens of percent. 



3.3. Generic Representation 

The squared matrix element Eq. (|T6|) involves a factor uo~ 2 which survives the nucleon 
phase-space integration and which, in fact, is a generic feature of any bremsstrahlung 
process (Raffelt 1996). Therefore, we will always write the lowest-order scattering kernel as 

S£\u) = T -^s{u/T) (19) 

where r CT is what we call the "spin fluctuation rate" while s(x) is a dimensionless, slowly 
varying function of order unity. The factorization between T a and s(x) is not unique. We 
take r CT such that s(0) = 1 when the neutrons are nondegenerate and when the pion mass 
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has been ignored in the matrix element, leading to (Raffelt and Seckel 1995) 



3tt 2 




T 2 



(20) 



rv 



m N 



In terms of the neutron Fermi momentum pp the "effective degeneracy parameter" is 



it varies relatively little throughout a SN core and thus is a convenient measure of the 
neutron density All modifications of Sjfi(u) by neutron degeneracy and the finite pion 
mass are included in the dimensionless function s(x) which will be determined in Sec. 4. 

For bremsstrahlung processes the singular behavior of S^(u) at uj = is of no concern 
because it is suppressed by a sufficiently high power of uj in all relevant neutrino phase-space 
integrations. However, we want a kernel which consistently describes bremsstrahlung and 
the "inelasticity" of scattering. Even in this case the singularity can be interpreted under 
a phase-space integral such that one obtains finite and meaningful results (Sawyer 1995; 
Raffelt, Seckel & Sigl 1996). Still, the singular behavior is an artifact of the lowest-order 
perturbative expansion. Multiple-scattering effects (or formally a resummation of the 
neutron propagator) render S a (oj) a well-behaved function everywhere. 

The interpretation of S a (uj) as a spin autocorrelation function suggests a 
phenomenological inclusion of multiple-scattering effects by the ansatz 



In the classical limit of "hard collisions" one has s(x) = 1 and r = IV, leading to the correct 
normalization. In the general case we choose T such that S a (u) is properly normalized. 

To summarize, our scattering kernel is determined by a two-step procedure. We first 
calculate the perturbative S^(lj) by what amounts to a bremsstrahlung calculation and 
from which we extract s(x). Next, we determine V in Eq. (^) such that the normalization 
Eq. ( |i~5|) is fulfilled. We stress that the exact shape of S a (uj) around uj = is not crucial; in 
that sense our Lorentzian ansatz is well-motivated but arbitrary. The important quality of 
our ansatz is that it gives the correct normalization and thus the correct neutrino scattering 
cross section. 
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3.4. Neutrino Processes 
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3. 4-1- Boltzmann Collision Integral 



We have defined the scattering kernel essentially as the squared matrix element of a 
neutrino- nucleon process with the nucleon degrees of freedom integrated out. In order to 
calculate quantities like a neutrino mean free path it remains to integrate over the neutrino 
phase space. To this end we begin with the Boltzmann collision integral for neutrinos 



where fi j2 are now the occupation numbers of neutrinos with momenta ki and k2, 
respectively, while fi t2 are those of antineutrinos with the corresponding momenta. The first 
two terms in square brackets represent pair emission and pair absorption by the medium, 
i.e. the bremsstrahlung process and its inverse, while the third and fourth term represent 
the gain and loss terms from neutrino scattering on the nuclear medium. An analogous 
equation obtains for antineutrinos, i.e. for fi. 

For a practical numerical implementation we mention that the dimension of S a is 
(energy) -1 so that an integral / du S a (u) is a dimensionless number. Further, in natural 
units we may write d 3 k 2 = dn 2 ^ 2 di^2 so that in the units usually employed in SN physics the 
integral expression in Eq. ( p3|) has the dimension MeV 2 . In the overall coefficient we write 
the baryon density as % = p/rn^ in terms of the mass density and the vacuum nucleon 
mass. Then the overall coefficient is numerically C\G 2 F nB = 3.79 x 10 4 s _1 MeV -2 pu, 
where p i4 = p/10 14 gcm~ 3 and where we have used Ca = —1.26/2. 



To illustrate the use of the scattering kernel and for later reference we turn to a few 
simple quantities related to bremsstrahlung. The inverse mean free path of a neutrino with 
energy u}\ against pair-absorption is 



fx 



coll 



C\G\n B J ^(3- cos B) 

x "(1 - /i)(l - h) S.i-uJr - u 2 ) - fJzSrfa + u 2 ) 

+ (1 - h)f 2 S^-Wi + u 2 ) -Ml- f 2 ) Srfa - UJ 2 )} , (23) 



3.4-2. Bremsstrahlung 



A 



-l 




/ 



d 3 k 2 



(3 - cos6>) f 2 S a (uj 1 + uj 2 ) 



(24) 



(2tt) 3 



with f 2 the occupation number of the antineutrinos. We use a Maxwell-Boltzmann 
distribution at temperature T for the neutrinos. With the representation Eq. (|19D for the 
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scattering kernel the spectral average is 



r^p M „„^nn — (Aj *) — — - n B T 2 ^ £ (25) 

7T Z(J7T 



with the "dimensionless mean free path" 

roc roc ™2™2 

£ = 5 / da* / dz 2 e-^e"* 3 7 ^-r s(xi + s 2 ). (26) 

Here, x% = uj\/T and x 2 = U2/T are the dimensionless neutrino and antineutrino energies 
respectively and we have ignored multiple-scattering effects. One of the energy integrations 
can be done analytically so that 

£ = -/ dxx 3 e~ x s(x). (27) 



6 

Taking s(x) = 1 leads to £ = 1; for realistic conditions it is in the range 0.2-0.5 (Sec. 4.5). 

Next we consider the energy-loss rate of a medium which is taken to be transparent to 
neutrinos so that one may ignore phase-space blocking effects, 

f d 3 k d 3 k 

Qnn^nnuu = C\ G F n B / TTrTTTTT^ ( 3 ~~ COS d ) ( W l + ^) S a (~UJi - L0 2 ) 

= ^ t r <t " /e " /T * (w) - <28) 

The production rate of neutrino pairs instead of the energy-loss rate is found by dropping 
the energy to = + uj 2 ) in this expression. 

In numerical SN codes it is sometimes useful to consider the differential energy 
production rate into the neutrino channel only, not counting the energy which goes into 
antineutrinos. Therefore, we must integrate over the v phase-space, ignoring for simplicity 
Pauli blocking effects, 



roc 

■u\ du{u-u x ) 2 e-" /T S a (u). (29) 



dw x 4tt 4 

The integral over this quantity is half of Qnn^nnvv of Eq. ( p8|) because we now measure 
only the energy carried away by the neutrino, ignoring that carried by the antineutrino. In 
proper units we have C\G f Ub = 1.13 x 10 31 ergs cm -3 s -1 MeV -6 pi±. 

In order to compare the efficiency of energy transfer of different processes it will 
turn out to be useful to consider the net transfer of energy between a neutron fluid at 
temperature T to a neutrino fluid at temperature T v . It is found to be 

AQ„„, = r dco^ (e~^ T - e-«/ T ») S». (30) 

407T 4 JO ' 
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With AT = T - T v and in the limit |AT| < T this is 



AQ 

AT 40vr 4 T 2 Jo 



3C 2 A G F n B T* r x 5 e~* 
- ^ Wo dx ^~ s ^ ( 31 ) 

where we have used the representation Eq. ( |I9|) for the lowest-order scattering kernel which 
may be used because the high power of u appearing under the integral renders the low-u; 
modification by multiple-scattering irrelevant. With s(x) = 1 the integral is 1. 



3.4-3. Neutrino Scattering 

When studying the bremsstrahlung process one is naturally led to the Boltzmann 
collision integral of Eq. fl23Q which would be incomplete without the scattering processes. 
For a neutrino of energy uii the scattering cross section differential with regard to the 
final-state energy 0J2 and the scattering angle is easily identified to be 

do-\ _ 3 - cos 9 C\G 2 F 2 S a {ux - u 2 ) ^ 

Numerically we have C\G 2 F = 2.10 x 10 -44 cm 2 MeV~ 2 . If the medium is dilute so that the 
spin-fluctuation rate is small, then S a {u) /2ir is strongly peaked around u = and actually 
must approach a 5-function which is normalized according to Eq. (|T5|). Integrating over 
energy and angles then gives us the usual elastic scattering cross section <Ti = {?>C\G 2 F / tt)uj\ 
which in a degenerate medium is reduced according to the right-hand side of Eq. ( PT5[ ) . 

Next we calculate the total axial-current cross section, averaged over a Maxwell- 
Boltzmann neutrino distribution. We find, in agreement with Raffelt & Seckel (1995), 

7T Jo livr 

In the dilute limit where S a (u) — > 2ttS(lj) the integral is 1. In general it is less than 
1 because the function S a (u) is normalized to 1 while the rest of the integrand falls 
monotonically with uj so that the integral has its maximum value when S a is narrowly 
peaked around uo = 0. Therefore, iViV interactions lead to a reduction of the average 
neutrino scattering cross section. This cross-section reduction is an unavoidable consequence 
of NN interactions and occurs on the same level of approximation as the bremsstrahlung 
emission or absorption of neutrino pairs. 
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In order to compare the efficiency of energy transfer of the "inelastic scattering process" 
with that of bremsstrahlung and other processes we finally consider the rate of energy 
transfer between a nucleon and a neutrino fluid at slightly different temperatures. We find 

AQ vnn ^ nnu 30 C 2 A G%n B T* r°° , 12 + 6x + x 2 



AT 



a Jo 20 6 ^ ^ 



where we have used Eq. ( |19"D for the lowest-order scattering kernel. With s(x) = 1 the 
integral is 1. 



4. Analytic Fitting Formula 

4.1. General Representation 

In order to represent the scattering kernel by analytic fitting formulae we write it in 
the form 

S ^ ) = T x 2 + { ig /2f s{x) (35) 

where x = to/T. The dimensionless quantities 7, g, and s(x) depend on the medium density 
and temperature. We characterize the former by the neutron effective degeneracy parameter 
introduced in Eq. fl2~l|) which may be written as 

^ = VU = 3-04 Pl4 T 10 . (36) 



The spin-fluctuation rate of Eq. ([20]) is numerically 

7 = jr = 1-63 ^/ 2 T 10 = 8.6 Pl4 Tr 1/2 . (37) 

Both g and s(x) are always of order unity so that simple estimates of bremsstrahlung- related 
quantities may be obtained by setting them both equal to 1. 

The matrix element Eq. (|1^) involves the pion mass as one more dimensionful 
parameter. It proves useful to express it in the form 

2 

TfL 

^^ = L94r - (38) 

With this notation we have that the dimensionless scattering kernel s(x) and the quantity 
g depend on the parameters y and 77*. 
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4.2. Nondegenerate Limit 



We will construct a general analytic representation from an interpolation between the 
nondegenerate (ND) and degenerate (D) limiting cases. The former was treated by Raffelt 
& Seckel (1995); in a symmetric form their result can be expressed as 



_ , \ s ND (x,y) + s ND (-x,y) 
SND{x,y) = 



e -x/2 _|_ e x/2 ^oo 



16 



roc 
J\x\ 



dt e 



_ t/2 3x 2 + Gty + W 
3 



2Vt 2 



x' 



log 



't + y + Vt 2 



x' 



x 2 + 2ty + y 2 t + y \t + y - ^t 2 - x 



(39) 



Evidently, this expression is even in x; the structure function is recovered by 
snd(^, y) = 2snd(^, y)/(l + e~ x ) and then obeys detailed balance. Its limiting 
behavior is 

'1 for x = and y — 0, 

SND(z,y) = S (47T/X) 1 / 2 for x > 1 and x > y, 

fy- 2 for a; = and y > 1, 

i.e. it falls off to zero for either large x or y. 



(40) 



In order to construct an analytical fit we first extract the detailed-balance behavior 
explicitly by 

for x > 0, 



s(x,y) = s(x,y) x 



for x < 0. 



(41) 



An analytical fit for the y = case (vanishing pion mass) is 



X 

4:71 



1+ (12 + -]a; 

7T 



-1/12N 



-1/2 



(42) 



This function reproduces the limiting behavior at x = and x ^> 1, it has the correct 
derivative 1/2 at x = 0, and it deviates from the correct result by no more than 2.5% 
anywhere. Likewise we have constructed 



-o(v) 



1 , V\ 
1+ 2 l0g 



y 



l 



,160 J 



2 \ V9' 



-9/4 



(43) 



It has the correct limiting behavior at y = and y ^> 1, it is vertical at y = in agreement 
with the full result, and its maximum error is below 4.5% anywhere. This function gives us 
a good idea of the suppression caused by including the pion mass. 

A simple estimate of the compound scattering kernel is snd(^,1/) — s y=0 (x) s x=0 (y). 
However, it underestimates the correct value by as much as a factor of a few when both 
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x and y are a few which is the most relevant regime in a SN core. To construct a more 
suitable approximation we expand the square bracket in Eq. ([39]) in a power series in t 
and keep only the lowest term which is proportional to t 3//2 . Now the exponential can be 
integrated so that 

27r 1 /V/ 2 (3a; 2 + 6xy + 5y 2 ^ 

valid for i>1 and t/> 1. With a little bit of tinkering one can supplement it such that it 
behaves reasonably for small x and y, 

27r V2 f x + 2- e-y/ 12 f /2 (x 2 + 2xy + \y 2 + l) 

W*, v) 7rl/2 + (7r i/B + x + , ) 4 • ( 45 ) 

This function is 1 for x = y = and has the correct limiting behavior for y = and x ^> 1, 
but it does not have the correct limiting behavior for x = and y ^> 1. The maximum error 
for any x is 14%, 7%, 5%, 9%, and 12% for y = 0.5, 1, 2, 3, and 5, respectively. The errors 
are even smaller for an integrated quantity like the average mean free path of Eq. (|2~T| ) where 
we find a deviation of —4%, —1%, 4%, 7% and 10% for the same y values. This precision is 
good enough in view of the relatively crude treatment of the nuclear matrix element. 



4.3. Degenerate Limit 

In the degenerate limit the phase space integration can be carried out analytically even 
with a nonzero pion mass. Ishizuka & Yoshimura (1990) found 

where 

*, s 1 5u f 2 \ u2 u2 /2 v / 2w 2 + 4\ 

t[u) = 1 arctan — M tt, r H ; arctan . (47) 

JK J 6 \uJ 3(u 2 + 4) 6V2u 2 + 4 \ u 2 ) y 1 

The function sn(x) fulfills detailed balance explicitly. However, it is pathological in that it 
scales as x 3 for x ^> 1 which implies that the full scattering kernel is not normalizable in 
the sense of Eq. (|T5|). Of course, the assumption of degenerate nucleons is never good for 
energy transfers which far exceed the nucleon Fermi energy. For u 3> Ep (x ^> Ep/T) the 
full s{x) must always approach snd(^)- 
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4.4. Interpolation 

An interpolation which is accurate to within roughly 30-40% for all values of x and y 
is provided by 

*=(— + -) ( 48 ) 

One can do better by using 

* = (%D fe) + sn (y) ) ~ lMy) F(x, !/, V*)[l + C(x, y, v*) G(x, y, ^)] , (49) 



where 



F(x,y,r}*) = 1 + 



[ 3 + ( x _l. 2 )2 + a: -4]( 1 + ?7 2)( 1 + !/ 4)> 
G(x,y,ri*) = 1 - 0.0044 x 11 



n ^ 

y v* 



C(x,y,r}*) 



0.8 + 0.06y 105 + 0.2' 
1.1 x 1A h(r]*) 30 



2.3 + h(r]*)x - 93 + 0.0001 x 1 - 2 30 + 0.005 x 2 - 8 ' 
p( y ) = 0.67 + 0.18 y - 4 , 

A <"«» - 2i5W- (50) 

This interpolation function reproduces the true value of the scattering kernel to within 
5-10% in the physically interesting parameter space, y ~ 2-6 and 77* ~ 2-6, and maximum 
deviations for any value of the parameters of less than 25%. As discussed in Sec. 3.3, the 
uncertainties from the nuclear matrix element are at least of that magnitude so that this 
fitting accuracy is entirely sufficient. 

In Table 2 we give s(x, y, 77*), both the numerically calculated values and for comparison 
the analytical fit. In Fig. 3 we plot s(x) for different values of y and 77* where the solid 
lines represent the numerical result and the short-dashed ones the fitting formula. The 
detailed-balance condition tells us that s(—x) = s(x)e~ x . Still, we plot s(x) somewhat 
across x = in order to illustrate that it is smooth at the origin. In both panels the 
upper curve corresponds to a vanishing pion mass (y = 0) and essentially nondegenerate 
conditions (77* = 0.31). 

The upper panel illustrates the effect of degeneracy which suppresses the scattering 
kernel at low x, but enhances it at large ones. All of these curves approach the nondegenerate 
case asymptotically in the limit x — > 00, but they do so from above. The lower panel 
illustrates the impact of including the pion mass in the matrix element. It suppresses the 
scattering kernel at low x by an amount which depends on y, and again the curves approach 
the y = case in the large x limit, however here they do so from below. 
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4.5. Multiple Scattering 

It remains to determine the function g(y,T]*) for the denominator of Eq. (^) which 
must be adjusted to meet the normalization condition Eq. (pi|), 

l-oo dx 7 

x , 2 + [ ig /2y s{x) = B - (51) 

For ND conditions and when 7CI the Lorentzian is essentially a 5 function and we have 
for the integral snd(0)/<7. Therefore, in the limit y = and 77* = 0, which also implies 
7 = 0, we have g = 1 because 7 was defined such that in this limit snd(0) = 1. 

We have determined the function g(y,r)*) using the normalization condition Eq. (|15|) 
and our analytic approximation formula Eq. ( |4"9"D for s(x). Furthermore we need to provide 
an analytic approximation for g(y,rj*) that can be implemented in numerical calculations. 
For very high values of y, we know that g oc y~ 2 and for y — > we know that g — > 0.5 in 
the degenerate limit and g — > 1 in the nondegenerate limit. A decent analytic fit that has 
the correct limiting behavior is given by 

9(y^V*) = ^ — Pl +2/io 52 

1 + a 3 y P2 + a 2 y Pl+ / 13.75 

where the coefficients are given by 

ai = — _ 1 2 ; - + (O.5 + 77/I5.6) 



a 2 



l + r/- 1 25y 2 + l v " '25y 2 + V 
0.63 + 0.04 r/ L45 



1 + 0.02 /7 2 - 5 ' 
a 3 = 1.2 exp(0.6r/ - 0.4/7 15 ) 
1.8 + 0.45/7 



Pi 



1 + 0.1577 1 - 5 ' 



0.05/7 / s 

p 2 = 2.3 — . 53 

y2 1 + 0.025/7 v ; 

In Fig. 4 we show g(y, 77*) as a function of y for different values of the degeneracy parameter 
77* (solid lines) as well as our analytic fit to g(y,r]*) (dashes). Also, in Fig. 5 we show the 
normalization -8(77*) which was defined in Eq. ( |T5D as a function of 77*. 

We may now calculate the average neutrino absorption rate in its dimensionless form 
Eq. (p7| ). In Table 3 we give values for £ based on our analytic approximation scheme, both 
with and without the inclusion of multiple-scattering effects. For realistic SN conditions 
where y ~ 2-6 and 77* ~ 2-6 we see that £ ~ 0.2-0.5. Multiple scattering is not a strong 
effect for the average absorption rate. 
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However, to calculate inelastic scattering it is unavoidable to include the multiple 
scattering effect. With our function g(y, r)*) the scattering kernel is self-consistent in that 
it allows for the simultaneous treatment of bremsstrahlung and inelastic scattering while 
reproducing the correct total scattering cross section. 

5. Conclusions 

We have studied the neutron-neutron bremsstrahlung process nn <-> nnuu as an opacity 
source for neutrinos and its impact on the neutrino spectra formation in supernovae. For 
and u T , bremsstrahlung is by far the most important number- changing reaction, far more 
important than e + e~ annihilation or the plasma process. As an energy-changing process 
it is roughly as strong as elastic scattering on electrons and positrons. We estimate that 
including bremsstrahlung will reduce the temperature by more than 1 MeV. This is only 
a differential effect in the sense that the and v e spectral temperatures become closer 
by this amount. Their absolute changes can not be determined by our simple procedure. 
The overall magnitude of our differential effect is in agreement with Suzuki's (1991, 1993) 
numerical simulations. 

If one includes the bremsstrahlung process it is inconsistent to ignore the inelasticity of 
vn scattering which may be pictured as the "crossed bremsstrahlung process" unn — > nnv. 
This process cannot change the neutrino number density, but as a source of spectral 
equilibration it is roughly a factor of 10 more important than bremsstrahlung and thus also 
far more important than ve scattering. Moreover, depending on details of the temperature 
and density profile of the SN core surface layers, the nucleon recoil in vN collisions is of 
roughly equal importance to the inelasticity (Janka et al. 1996). Bremsstrahlung, inelastic 
scattering, and recoils near the neutrino spheres is likely to make the spectra and fluxes of 
the different neutrino species far more similar than had been thought previously. This would 
modify all phenomena related to neutrino flavor oscillations which rely on the difference in 
the spectral temperatures to be effective. 

The main purpose of our paper was to provide an explicit form for the scattering kernel 
S a {ui) that will allow one to study the impact of bremsstrahlung and inelastic scattering 
consistently in a numerical simulation. 

Our derivation is fundamentally perturbative; it is useful to study the spectra 
formation in the surface layers of a proto neutron star. We believe our scattering kernel 
to be appropriate for densities roughly below 10 14 gcm~ 3 . At higher densities our result 
predicts a huge suppression of the total uN scattering cross section which opens the thorny 
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issue of its actual magnitude. This problem has not been solved so that any SN simulation 
has to rely on an arbitrary prescription for the neutrino diffusion coefficients in the deep 
interior. The SN 1987A signal (Keil, Janka & Raffelt 1995) as well as the /-sum rule for 
the spin-density structure function (Sigl 1996) indicate that the cross-section suppression 
is not as large as predicted by our scattering kernel, but how large the cross sections truly 
are is not known. Recent attempts to include the correct nucleon dispersion relation as well 
as hyperons (Reddy & Prakash 1997) predict modifications of the "standard cross sections" 
by factors of order unity, but no attempt was made to address the modifications caused by 
the inevitable nucleon or hyperon spin-spin interactions. 

Another shortcoming of our scattering kernel is that it does not include recoil effects. 
In some regions of a SN core it may be even more effective than NN interactions at 
modifying the neutrino energy in a collision. A scattering kernel which includes recoils and 
NN interactions simultaneously requires avoiding the long-wavelength approximation, a 
rather complicated task that would imply another variable (the scattering angle) explicitly 
in the scattering kernel. 

Our estimates clearly show that bremsstrahlung and its related processes cannot be 
ignored for the formation of the neutrino spectra in a SN core. However, to determine 
the real magnitude of the modification one needs to perform a self-consistent numerical 
analysis which includes feedback effects on the medium. If such a study reveals that the 
modifications are as severe as suggested by our estimates one may be motivated to derive a 
more complete perturbative scattering kernel. One could study the proton-neutron process 
on the basis of a realistic nucleon-nucleon potential in the spirit of Sigl's (1997) recent 
work, thereby avoiding the pathologies inherent in an OPE Born treatment of np scattering. 
Moreover, it would be cumbersome but straightforward to include recoil effects and NN 
interactions simultaneously in the scattering kernel. Unfortunately, it remains far from 
obvious how to proceed extending the scattering kernel into the high-density regime. 
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A. Nucleon Phase-Space Integration 

We give some details of solving the nucleon phase-space integral Eq. From the 

original 12- dimensional integral one can integrate out 8 dimensions analytically so that 
in the end only a four dimensional integral remains to be done numerically. First, the 
three-dimensional momentum delta function is used to integrate out d 3 p4. Because the 
nucleons are nonrelativistic their energy is E^ = pf/2m^ so that 

E 1+ E,-E,-E t + u= -**- 2pi - p 't 2pi - p > + 2p '- p > + ». (Al) 

2m N 

Because of the assumed isotropy of the medium the pi momentum may be chosen in the 
^-direction, i.e. we have trivially / <i 3 pi = Am J dp\ with p\ = |pi|. Next we use polar 
coordinates with a and (3 the polar and azimuthal angles of P2 relative to pi and 9 and 4> 
those of p 3 so that 

<i 3 p2 = pi dp2 d cos a d(5 (A2) 
d 3 ps = pi dps d cos 9 dcf). (A3) 

Because of the medium's isotropy one of the azimuthal integrations is trivial. We chose 
the dcj) integration so that three nontrivial angular integrations remain. In terms of these 
remaining angular variables we have 

Pi P2 = P1P2 cos a (A4) 
Pi • P3 = PiP?, cos 9 (A5) 
P2 • P3 = cos a cos 9 + sin a sin 9 cos f3. (A6) 



The integration over d(3 is carried out using the 5-function where we write 
f{(3) = E 1 + E 2 - E 3 - Ei + u so that 



r27r p2tt I 

j Q dp6(f(P)) = Y,J 8{p-fr) 



dfW) 




d(5 





(A7) 



where (3i with i — 1, 2 are the two roots of f{0) = 0, one in the interval [0, tt] and one in 
the interval [71*, 2ir]. Because Eq. (|A^7|) is symmetric in j3 it is enough to use fii which is in 
the interval [0, 7r] and multiply by 2. Thus we end up with 



2tt 



d(58{f{(5)) 







d/3 





e 



df(p) 



dp 



(A? 



/3=/3i 
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The derivative may be written in the form 

df((3) 



dp 



= Vaz 2 + bz + c, (A9) 

13=13! 



where z = cos a and 



a = pl(- P l-pl + 2 PlP3 cose), 

b = 2uom N pip 2 - 2pip 2 P3 - 2ujm N p 2 p 3 cos 9 + 2p\p 2 p 3 cos 9 
+ 2p 2 p\ cos 9 — 2p\p 2 p\ cos 2 9, 

O O O /I „ 

c = a; mjy + 2umNp 3 + p 2 p 3 — p 3 — 2uomNpip 3 cos 9 

+ 2pip| cos 6* - p\p\ cos 2 6> - p\p\ cos 2 6>. (A10) 

The integration of the 5 function also fixes E 4 = E 1 + E 2 — E 3 + uj which is used to 
determine the Pauli-blocking factor of particle 4 in the phase-space integral. 

In order to perform the integration over dz = d cos a analytically we note that partial 
integrations allow us to reduce the integrand to a sum of two terms. One of them is 
independent of z while the other is of the form l/(qz + r) where q and r are expressions 
which do not depend on z. We find explicitly 



Q(az 2 + bz + c) = —=Q{b 2 -Aac), 



x a/ az 2 + bz + c v —a 



f 

J —i 



/ 

- Q(az 2 + bz + c) = = \. = \ All 

^/az 2 + bz + c qz + r ^-ar 2 q- 2 + brq~ l - c 



The step function in the integrand singles out the physically relevant range where 
-1 < z < +1. 

After these analytic manipulations we are left with a four dimensional numerical 
integral over dp\dp 2 dp 3 d cos 9 where the integration region is determined by the step 
function 0(6 2 — 4ac). 
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Table 1. Neutrino Energy Spheres 



Flavor: 

/brcms 





1 


5 


v T 
10 


20 





1 


5 


10 


20 


p sph [10 14 g/cm 3 ] 
T sph [MeV] 
r~ h [MeV] 

tot 
T sph 


1.20 
8.78 
7.39 
16.3 


0.99 
7.50 
6.32 
7.11 


0.85 
6.67 
5.63 
4.42 


0.78 
6.31 
5.33 
3.49 


0.71 
5.99 
5.06 
2.83 


0.82 
6.52 
5.49 
3.83 


0.77 
6.25 
5.25 
2.78 


0.70 
5.94 
5.02 
2.18 


0.66 
5.77 
4.87 
1.93 


0.62 
5.60 
4.73 
1.76 
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R [km] 

Fig. 1. — Physical parameters of the proto neutron star model S2BH_0 of Keil, Janka & 
Raffelt (1995) which represents a SN core 1 s after bounce. For the neutrons we show the 
effective degeneracy parameter 77* as defined in Eq. (^1|), taking the vacuum nucleon mass 
for win- 
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V_J I I I I I 

10 20 30 40 50 
E„ [MeV] 

Fig. 2. — Temperature of the medium (SN model of Fig. 1) at last energy exchange for 
different processes: Scattering on e ± (solid line), inverse bremsstrahlung vvnn — > nn (short 
dashes), and /5-absorption u e p — > ne (long dashes). 
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5 10 15 
x 

Fig. 3. — Dimensionless scattering kernel s(x, y, rj*). Upper Panel: Pion mass ignored (y = 0) 
with r]* = 0.31, 1.01, and 2.42 from top to bottom, corresponding to r] — —2, 0, and 2, 
respectively. The solid curves are numerical results whereas the short-dashed ones are the 
analytical approximation Eq. fl49"l) . Lower Panel: Nondegenerate case (rj = —2, 77* = 0.31) 
with y — 0, 2, and 4 from top to bottom. 
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Fig. 4. — Function g(y,r}*) as denned in Sees. 4.1 and 4.5 to take multiple-scattering effects 
into account. From bottom to top the curves are for 77* = 0.31, 1.01, 2.42, 4.22, 6.14, and 
8.11 respectively. The numerically calculated g(y,T]*) are shown as full lines, whereas the 
analytic fit from Eq. Q52"D is represented by the dotted lines. 
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Fig. 5. — The normalization function B(i]*) defined in Eq. (|15[) as a function of the nucleon 
degeneracy parameter 77*. 
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Table 2. s(x,y,r)*) in the form "Numerical Value/ Analytical Approximation". 



y = 





2 4 


6 8 


j] = (77, = 1.01) 



x = 0.701/0.735 0.267/0.277 0.166/0.171 0.117/0.117 0.088/0.086 

1 0.861/0.799 0.389/0.412 0.248/0.272 0.177/0.193 0.134/0.145 

2 0.926/0.916 0.482/0.508 0.321/0.352 0.233/0.259 0.179/0.198 

3 0.959/0.955 0.551/0.579 0.381/0.416 0.284/0.314 0.222/0.246 

4 0.973/0.977 0.600/0.633 0.429/0.466 0.327/0.360 0.259/0.287 

5 0.976/0.991 0.636/0.673 0.467/0.506 0.362/0.398 0.291/0.322 

6 0.971/0.998 0.660/0.702 0.496/0.538 0.391/0.429 0.319/0.351 

7 0.959/1.000 0.675/0.724 0.518/0.562 0.414/0.454 0.341/0.375 

8 0.945/0.998 0.685/0.738 0.534/0.581 0.433/0.474 0.360/0.395 

9 0.928/0.993 0.690/0.748 0.546/0.595 0.448/0.489 0.376/0.411 
10 0.911/0.985 0.691/0.753 0.554/0.605 0.459/0.501 0.389/0.424 



77 = 2 (77, = 2.42) 



= 0.349/0.360 0.165/0.172 0.111/0.115 0.082/0.084 0.064/0.064 

1 0.485/0.443 0.252/0.270 0.172/0.190 0.128/0.141 0.101/0.110 

2 0.600/0.567 0.342/0.368 0.240/0.267 0.182/0.203 0.144/0.160 

3 0.705/0.676 0.431/0.463 0.310/0.343 0.238/0.266 0.191/0.212 

4 0.794/0.770 0.513/0.550 0.378/0.415 0.295/0.326 0.238/0.263 

5 0.867/0.847 0.585/0.624 0.440/0.478 0.349/0.381 0.285/0.311 

6 0.924/0.909 0.646/0.685 0.495/0.532 0.397/0.429 0.328/0.353 

7 0.964/0.957 0.694/0.735 0.541/0.577 0.439/0.469 0.366/0.390 

8 0.993/0.992 0.733/0.773 0.580/0.613 0.476/0.502 0.400/0.421 

9 1.011/1.016 0.763/0.801 0.611/0.641 0.506/0.529 0.428/0.446 
10 1.020/1.032 0.784/0.821 0.636/0.663 0.531/0.550 0.453/0.466 



77 = 4 (77* = 4.22) 



= 0.153/0.154 0.086/0.087 0.063/0.064 0.049/0.050 0.040/0.040 

1 0.229/0.219 0.136/0.145 0.100/0.109 0.078/0.086 0.064/0.070 

2 0.313/0.312 0.197/0.218 0.147/0.167 0.117/0.133 0.096/0.109 

3 0.408/0.415 0.269/0.304 0.204/0.235 0.163/0.189 0.135/0.156 

4 0.509/0.522 0.349/0.396 0.269/0.310 0.217/0.252 0.181/0.209 

5 0.613/0.625 0.434/0.489 0.339/0.387 0.276/0.316 0.231/0.263 

6 0.712/0.720 0.518/0.577 0.409/0.460 0.337/0.378 0.284/0.317 

7 0.803/0.802 0.598/0.654 0.478/0.525 0.397/0.434 0.336/0.366 

8 0.884/0.871 0.671/0.720 0.543/0.582 0.454/0.484 0.387/0.410 

9 0.953/0.927 0.736/0.773 0.602/0.629 0.506/0.525 0.435/0.447 
10 1.011/0.970 0.793/0.814 0.654/0.667 0.554/0.559 0.478/0.477 
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Table 2 — Continued 



V = 





2 4 


6 8 


r) = 6 (77* = 6.14) 



= 0.074/0.074 0.046/0.046 0.036/0.036 0.029/0.029 0.025/0.025 

1 0.114/0.116 0.074/0.080 0.057/0.063 0.047/0.052 0.040/0.044 

2 0.164/0.177 0.111/0.127 0.087/0.101 0.072/0.084 0.061/0.071 

3 0.227/0.253 0.159/0.189 0.126/0.152 0.104/0.126 0.088/0.107 

4 0.300/0.340 0.217/0.265 0.173/0.214 0.144/0.179 0.123/0.152 

5 0.385/0.435 0.284/0.350 0.229/0.284 0.192/0.238 0.165/0.203 

6 0.476/0.530 0.359/0.438 0.293/0.359 0.247/0.301 0.212/0.258 

7 0.572/0.621 0.440/0.526 0.361/0.432 0.306/0.364 0.264/0.313 

8 0.669/0.705 0.522/0.606 0.432/0.501 0.368/0.424 0.319/0.364 

9 0.763/0.779 0.603/0.678 0.503/0.563 0.431/0.477 0.375/0.411 
10 0.851/0.841 0.681/0.738 0.572/0.615 0.492/0.523 0.431/0.452 



7/ = 8 (ry* = 8.11) 



= 0.043/0.041 0.029/0.027 0.023/0.022 0.020/0.018 0.017/0.016 

1 0.068/0.068 0.046/0.048 0.037/0.039 0.031/0.033 0.027/0.029 

2 0.100/0.108 0.071/0.079 0.057/0.065 0.048/0.055 0.042/0.048 

3 0.142/0.161 0.103/0.122 0.084/0.101 0.072/0.086 0.062/0.074 

4 0.194/0.227 0.145/0.178 0.119/0.147 0.102/0.126 0.088/0.109 

5 0.257/0.302 0.196/0.246 0.163/0.204 0.139/0.175 0.121/0.152 

6 0.331/0.385 0.257/0.322 0.214/0.269 0.184/0.230 0.161/0.201 

7 0.413/0.469 0.325/0.403 0.273/0.338 0.236/0.290 0.207/0.253 

8 0.503/0.552 0.401/0.484 0.339/0.408 0.294/0.351 0.259/0.306 

9 0.597/0.630 0.482/0.561 0.409/0.474 0.356/0.409 0.315/0.357 
10 0.693/0.700 0.566/0.631 0.483/0.535 0.422/0.462 0.374/0.404 



Table 3. Dimensionless mean free path against pair-absorption £ as defined in Eq. 
Values are given as "without /with" multiple scattering. 



77* = 1 2 4 6 8 

y=l 0.750/0.712 0.666/0.488 0.471/0.102 0.321/0.017 0.222/0.005 

2 0.609/0.604 0.557/0.524 0.408/0.237 0.284/0.052 0.199/0.013 

4 0.447/0.447 0.417/0.415 0.320/0.297 0.230/0.157 0.166/0.063 

6 0.345/0.345 0.327/0.327 0.259/0.256 0.192/0.179 0.141/0.112 

8 0.276/0.276 0.264/0.264 0.215/0.215 0.164/0.161 0.123/0.116 



